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We study the renormalization group flow in a class of scalar-tensor theories involving at most two 
derivatives of the fields. We show in general that minimal coupling is self consistent, in the sense 
that when the scalar self couplings are switched off, their beta functions also vanish. Complete, 
explicit beta functions that could be applied to a variety of cosmological models are given in a five 
parameter truncation of the theory in d = 4. In any dimension d > 2 we find that the flow has only 
a "Gaussian Matter" fixed point, where all scalar self interactions vanish but Newton's constant 
and the cosmological constant are nontrivial. The properties of these fixed points can be studied 
algebraically to some extent. In d = 3 we also find a gravitationally dressed version of the Wilson- 
Fisher fixed point, but it seems to have unphysical properties. These findings are in accordance 
with the hypothesis that these theories are asymptotically safe. 



I. INTRODUCTION 

Fundamental scalar fields have not yet been observed, but they play a crucial role in the standard model and in grand 
unified theories, as the order parameters whose VEV is used to distinguish between otherwise undifferentiated gauge 
interactions. Whether such scalar order parameters are elementary fields, as in the standard model, or composites, 
as in technicolor theories, is still an open question. Known examples of the Higgs phenomenon (superconductivity, 
the chiral condensate in QCD) point to the latter possibility, but even if this was the case it might still be possible to 
use scalar theory as an effective description (a la Landau-Ginzburg) at sufficiently low energy. 

Scalar fields also play an important role in theories of gravity. Due to their simplicity they are very often used as 
models for matter. Also, because of the ease by which one can generate a nontrivial VEV, with an energy momentum 
tensor that resembles a cosmological constant, a scalar field is the most popular option as a driver of inflation. 
Furthermore, scalar fields easily mingle with the metric: by means of Weyl transformations it is possible to rewrite 
the dynamics in different ways sometimes leading to new insight or to simplifications. Theories of gravity in which 
a scalar is present are often called scalar-tensor theories. In this paper we will discuss the quantum properties of a 
class of theories of this type. 

The original motivation for this work comes from the progress that has been made in recent years to- 
wards understanding the UV behaviour of gravity. It seems that pure gravity possesses a Fixed Point (FP) 
with the right properties to make it as ymp totically safe, or in other words nonperturbatively renormalizable 
II II H, U 0, MM El El El El El El El, E3, El, El HI (see also for reviews). Let us assume for a 
moment that this ambitious goal can be achieved, and that pure gravity can be shown to be asymptotically safe. Still, 
from the point of view of phenomenology, we could not be satisfied because the real world contains also dozens of 
matter fields that interact in other ways than gravitationally, and a their presence affects also the quantum properties 
of the gravitational field, as is known since long [2l|. Indeed, in a first investigation along these lines, it was shown 
in [22| that the presence of minimally coupled (i.e. non self interacting) matter fields shifts the position of the gravi- 
tational FP and the corresponding critical exponents. In some cases the FP ceases to exist, so it was suggested that 
this could be used to put bounds on the number of matter fields of each spin. More generally the asymptotic safety 
program requires that the fully interacting theory of gravity and matter has a FP with the right properties. Given the 
bewildering number of possibilities, in the search for such a theory one needs some guid ing principle. One possibility 
that naturally suggests itself is that all matter self-interactions are asymptotically free [33]. Then, asymptotic safety 
requires the existence of a FP where the matter couplings approach zero in the UV, while the gravitational sector 
remains interacting. We will call such a FP a "Gaussian Matter FP" or GMFP. Following a time honored tradition, as 
a first step in this direction, scalar self interactions have been studied in [341. [3oT|. Here we pursue that study further. 
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The tool that we use is the Wetterich equation, an exact renormalization group (RG) flow equation for a type 
of Wilsonian effective action 1^, called the "effective average action". This functional, depending on an external 
energy scale k, can be formally defined by introducing an IR suppression in the functional integral for the modes with 
momenta lower than k. This amounts to modifying the propagator of all fields, leaving the interactions untouched. 
Then one can obtain a simple functional RG equation (FRGE) for the dependence of Tk on k (25l. l26l l27l. l28l. l29l. |30| . 
Insofar as the effective average action contains information about all the couplings in the theory, this functional 
RG equation contains all the beta functions of the theory. In certain approximations one can use this equation to 
reproduce the one loop beta functions, but in principle the information one can extract from it is nonperturbative, in 
the sense that is does not depend on the couplings being small. 

The most common way of approximating the FRGE is to do derivative expansion of effective average action and 
truncate it at some order. In the case of a scalar theory the lowest order of this expansion is the local potential 
approximation (LPA), where one retains a standard kinetic term plus a generic potential [25|, [26|, |3j ■ In the case 
of pure gravity, the derivative expansion involves operators that are powers of curvatures and derivatives thereof. 
This has been studied systematically up to terms with four derivatives [ll|, [HI, and for a limited class of operators 
(namely powers of the scalar curvature) up to sixteen derivatives of the metric[13l [T3 | . In the case of scalar tensor 
theories of gravity, one will have to expand both in derivatives of the metric and of the scalar field. 

In this paper we study (Euclidean) effective average actions of the form 



This can be seen as a generalization of the LPA, where one also includes terms with two derivatives of the metric. 

In [35[ it was shown that in d = 4, and assuming that V and F are polynomials in </> 2 , this theory admits a GMFP 
where only the lowest (</> 2 -indcpcndent) coefficient in V and F are nonzero. In this paper we extend and generalize 
this result in various ways. First of all, using so called "optimized" cutoff types [HI it is possible to write the beta 
functions in closed form, whereas in 13511 they could only be studied numerically. This makes the subsequent analyses 
much more transparent. Unlike in [35j, we will not assume from the outset that V and F are polynomials. Then, 
using the optimized cutoff it is possible to write explicit beta functionals for V and F . Exploiting general properties 
of these functionals we will be able to prove certain properties of the linearized flow in the neighborhood of the FP 
which had only been numerically observed previously. Namely we show that the matrix describing the linearized flow 
only has nonzero entries on the diagonal and on three lines next to it, and furthermore it has a block structure such 
that knowledge of the first two 2x2 blocks determines the whole matrix. 

The discussion in this paper is also more general than that of [35| in two ways: we keep the dimension of spacetime 
arbitrary and we allow for a more general gauge fixing, depending on two arbitrary parameters. Keeping the dimen- 
sion general is useful in view of possible applications to popular "large extra dimensions" theories Q, and also to 
higher dimensional dilatation-symmetric models which lead to vanishing cosmological constant in four dimensions [24| . 
Furthermore, with the closed form beta functions we can also perform a better search for other FP's where the scalar 
interactions are not all turned off. In [35| a numerical search was conducted on a grid of points in the neighborhood 
of the Gaussian FP, and no nontrivial FPs were found. Here, having the closed form of the beta functions, we can 
look for FPs by different methods. In d = 4, 5 and 6, some such points are found, but they appear to be spurious. 
On the other hand in d = 3 there is a FP which seems to be a genuine generalization of the Wilson-Fischer FP [Z(| , 
but it has unphysical properties. 

In a companion paper [4l| we will extend the discussion to a more general class of effective actions, 



where L is in general a function of the curvature scalar and of the scalar field. 

Even though the original motivation of our work was to study the UV properties of the theory, it is important to 
stress that the beta functions that we obtain are completely general: they hold for any energy range. Depending on 
the ratio between the parameters of the theory (the cosmological constant, Newton's constant, the scalar mass and 
all the dimensionful higher couplings) different terms in the beta functions will come to dominate. However, this is 
something that has not been put in a priori. Thus the beta functions can be used also to study IR or mesoscopic 
problems, provided the system can be accurately modelled by a scalar tensor theory. 

There is clearly much scope for applications to cosmology. Early work in this direction has been done in [23[ , using 
the beta functions of pure gravity. Along a different line, given the role played by scalar fields in inflation it seems 
likely that the RG running of couplings could have significant effects. This seems particularly true of recent attempts 
to use the standard model Higgs field as an inflaton, which use a special case of the action (QJ 
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with a large value for the nonminimal coupling £ [37j ]. The beta functions given in Appendix A. contain the full 
dependence on A, G, m 2 , A and £, including threshold effects and a resummation of infinitely many perturbative 
contributions. 

One can imagine also applications in the IR, for example along the lines of [24j , We mention that the appearance 
of a scalar field in the low energy description of gravity has been also stressed in (3|| . For a FRGE-based approach 
to that issue see also (39J. 

This paper is organized as follows. In section 2 we will derive the "beta functionals" for V and F. In section 3 we 
discuss the general properties of the GMFP, in any gauge and dimension. In section 4 we discuss numerical solutions 
for the GMFP. In section 5 we will discuss other FP's with nontrivial potentials and we conclude in section VI with 
some additional remarks. Apendix A contains some lengthy formulae for the beta functions of five couplings in four 
dimensions. 



II. THE BETA FUNCTIONS 



In this paper we will obtain beta functionals for the functions V and F defined in (JTJ) - To achieve this, we use 
Wetterich's functional renormalization group equation (FRGE) [25| 



d t T k = -STr 



5$ <5$ 



n k d t n k 



(3) 



where $ are all the fields present in the theory and STr is the generalized functional trace including a minus sign for 
fermionic variables and a factor 2 for complex variables, and lZ k is a suitable tensorial cutoff. 



A. Second variations 



In order to evaluate the r.h.s. of (J3]) we start from the second functional derivatives of the functional ([2]). These 
can be obtained by expanding the action to second order in the quantum fields around classical backgrounds: g pv = 
9p,u + hjxi, and (f> = <fi + 5(f), where <f> is constant. The gauge fixing action is given by 



(4) 



■XT'- 

and S g h is the corresponding ghost action given by 

S G h = - I d d x^C» 



£±Vft 



c n 



(5) 



These terms are already quadratic in the quantum fields. The second variation of eq. ([TJ) is, 



,(2) _ 



= \j d d xjg 



1.2 1 



-h 2 - iv/>"" j {v(f) - F^ 2 ) R) + F{4> 2 ) h + -F^ 2 ) ft □ ft 



-F^ 2 ) h^V^J p h p v - F(0 2 ) h a " h^ a iV - F(0 2 ) R pltm W° - -F{4> 2 )h^ □ V - F(0 2 ) hVpVvhf"' 



d d x^ 



ft ■ <f> (V'(<j) 2 ) - F'((f> 2 ) R) 8<j> + 2 cj>F'((j) 2 ) ft"" R^ 5<j> - 2 (f)F'((f> 2 ) 5(j> {V^ v h^ v - Uh) 



+ - j d d x^60 



- a + 2 V [<p 2 ) + 4 <i> 2 V " {(j) 2 ) - R (2 F' {(j) 2 ) + 4 cj) 2 F" [cp 2 )) 



Sgf + S, 



(6) 



Since we will never have to deal with the original metric g pv and scalar field <fi, in order to simplify the notation, in 
the preceding formula and everywhere else from now on we will remove the bars from the backgrounds. As explained 
in detail in Q, the functional that obeys the FRGE Q depends separately on the background field g pv and on a 
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"classical field" {g c i)av = g^v + {h c \) where {h c {) is Legcndre conjugate to the sources that couple linearly to 
h^. The same applies to the scalar field. In this paper, like in most of the literature on the subject, we will restrict 
ourselves to studying the effective average action in the case when (g c i) M „ = g^v and <j> c \ = <fi. From now on the 
notation g^ v and 4> will be used to denote equivalently the "classical fields" or the background fields. For a discussion 
of the effective average action of pure gravity in the more general case when (g c i) Mt/ ^ g^v we refer to [42j . 



B. Decomposition 



In order to partially diagonalize the kinetic operator, we use the decomposition of into irreducible components 



1 



1 



-flU 



(7) 



where hT v is the (spin 2) transverse and traceless tensor, £ M is the (spin 1) transverse vector component, a and h 
are (spin 0) scalars. In some cases this decomposition allows an exact inversion of the propagator. This happens 
for example in the case of maximally symmetric background metric. Thus with that in mind we will work on a 
c?-dimensional sphere. This change of variables in the functional integral gives rise to Jacobian determinants, which 

however can be absorbed by further field re-definitions ^ = y— □ — -f ^ and a = \/ —U^ —U — -^jo 0, @, H3]. 

Then the inverse propagators for various components of the field are easily read from the second variation of the 
effective action. Thus for the spin-2 component /i^„ we get the following inverse propagator: 



3d 



R 



d(d-l) 

For the spin-1 component £ we have the following inverse propagator: 



d-2 



F(4?)R. 



(8) 



(9) 



The two spin-0 components of the metric, a and h, mix with the fluctuation of <f> resulting in an inverse propagator 
given by a symmetric 3x3 matrix S, with the following entries: 



1 





1 










) 






Q 





























1, 



F{tf) (-□) - -V(tf) + 



2d 



F{<t> 2 )R 



F( 



□ -□ 



R 



d-1 




-□ -□- 



R 



2F{ ( f 2 ){-n) + V{cj?) 



1 

d-4 



F{cj> 2 ) R 



)(-a) + <pv'(<j ) 2 )- 0F'(0 2 )i? 

) + 4<j) 2 V"((j) 2 ) - (2F'(0 2 ) + 4(t> 2 F"((/) 2 ))R . (10) 

In order to diagonalize the kinetic operator occuring in the ghost action eq. ([5]) , we perform a similar decomposition 
of the ghost field into transverse and longitudinal parts in the following manner: 



c = c^ T + v^c , c„ = cj + 

where C^ T and Cj satisfy the following constraints, 

V M C" T = , v^cj = 0. 

Again this decomposition would give rise to a non trivial Jacobian in the path-integral, which is cancelled by the 
further redefinition C = \J — DC. For spin-1 component of the ghost, the inverse propagator is 



(11) 
(12) 



□ 



R 



(13) 
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while for spin-0 component we have the following inverse propagator 

<->-¥■ 

Now we have to specify the cutoff TZk occuring in FRGE eq. ([3]). We define TZk by the rule that T^ 2 ' + TZk has the 

(2) 

same form as T k ' except for the replacement of — □ by P/-(— □), where Pk{z) = z + Rk{z). Rk{z) is a profile function 

which tends to fc 2 for z — > and it approaches zero rapidly for z > k 2 . The quantity 1^ + ^fc i s the "modified 
inverse propagator". This procedure applies both to the bosonic degrees of freedom and to the ghosts. The cutoff 
TZk occuring in the FRGE depends on fc not only through the profile function Rk(z), but also through k dependent 
couplings present in the function F{<p 2 ) and F'(<f> 2 ). Thus the derivative fc^r = ^ acts not only on the profile function 
Rk(z), but also on the fc-dependent couplings present in F(cf> 2 ) and F'((j> 2 ). When this is neglected one recovers the 
one loop results. The presence of the beta functions on the RHS of the FRGE, produces a coupled system of linear 
equations, which has to be solved algebraically to yield the beta functions. 



C. The /3-functionals in d = 4 

To read off the beta functions we have to compare the r.h.s. of the FRGE with the i-derivative of eq. (Q}, namely 

d t T\g,<j>} = J d d x^g(d t V(cf> 2 ,R)-d t FR) . (15) 

(the gauge fixing and the kinetic term are not allowed to run in our approximations). Since the background R and 
(j> are constant, the space-time integral produces just a volume factor, which eventually cancels with the same factor 
appearing on the RHS of the FRGE. Thus the running of V and F can be calculated using, 



d t V{<i> 2 )= -^dtTk 



ftW) = -7777 



1 d{d t Yk) 



Vol OR 



(16) 



R=0 



where Vol is the space-time volume. In order to exhibit the explicit form of these beta functionals we go to d = 4, 
where Vol — 38 ^ , and set a = and /? = 1 (De-Donder gauge). Furthermore, we choose the optimized cutoff 
Rk{ z ) = {k 2 — z)9(k 2 — z) [36[, which allows to perform the integrations occuring in FRGE trace in closed form (see 
appendix A in [lJl). From the FRGE we then get 



d t V 



d,F 



1927T 2 1 

k 2 

2304tt 2 



150 



30^ 6(fc 2 * + 24 2 fc 2 F'*' + fc 2 FS 1 ) 

~gT + A 

\20k 2 F(ik 2 F -V) 24 



4 5 k 2 k 2 Tii 
F + + A 



240 2 fc 2 *' , ) 
d t F + — v — d t F' ) (17) 



* 2 



^ (24 cj) 2 k 2 F' + k 2 * + k 2 FEi) - 
A A 



A 



4 2 (6fc 4 F ,2 + V)A 



-4 <t> 2 f (7 k F' - V) (Si — k ) + 4 2 Si (7 k F' - V) (2 W V' - V *') + 2 fc <£ 2 S 2 + 48 fc 4 F' 2 * S 2 
9 t F 



-24 fc 4 ^ 2 Eg 



F 



30 



iofc 2 F(7* + 4y ) + 6 / 2FSlA + 4 ^ WA _ 24fc4j ^, 2s2 

w A \ 



4</> 2 k 2 F*'Si(7fc 2 F' - y') 



, 24 fc 2 



(fc 2 F' + 5 V) A - 12 fc 2 * E 2 - 2 (7 fc 2 F' - V) * Si 



(18) 



where we have defined the shorthands: 

* = fc 2 F-^; Si = k 2 + 2V' + A(f) 2 V" 



2 F' + 4 6 2 F" 



2 ,Tr'2 



*s/ 



Let us define the dimensionless fields = fc 2 (f), R = k 2 R and the dimensionless functions V^0 2 ) = fc d V^(0 2 ) 
and F(cf> 2 ) = k 2 ~ d F((j> 2 ). The beta functionals of the dimensionless and dimensionful functions are related as follows: 



[d t v)[4> 2 ] = -dv($ 2 ) + (d-2)4> 2 v'($ 2 ) + k- d (d t v)& 2 } , 

(d t F)[4> 2 } = -{d-2)F^ 2 ) + (d-2)^ 2 F'^ 2 ) + k-^(d t F)[cp 2 ] 



(19) 
(20) 
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Some comments are in order. From the expressions of dtV and dtF given in eq. (jTTJ) and (fTB)) respectively, we note 
that where ever there is occurence of (ft 2 , it occurs in combinations like 4> 2 V'V , (ft 2 V F' , (ft 2 F'F' , (ft 2 V'dtF' , (ft 2 F'd t F', 
<ft 2 V" and <ft 2 F" . Occurence of such combinations are crucial, as they help us (as is demostrated in |41() in proving 
that minimal coupling is self consitent. Because of the occurrence of dtF' in the r.h.s. of (fTT|) . the system of equations 
cannot be solved algebraically for dtF. It may be possible to solve it as a differential equation, but here we shall not 
pursue this. Rather, we observe that if F and V are assumed to be finite polynomials in (ft 2 of the form 

a b 

f(0 2 ) = ]Tw 2ri ; F(4> 2 ) = J2&n4> 2n . (21) 

n=0 n=0 

with finite a and b, then dtF' is also a finite polynomial in the beta functions and it becomes possible to solve for the 
beta functions algebraically. As an explicit example, in the appendix we give these equations in the de-Donder gauge 
(a = and 13=1) in d = 4 with five couplings truncation (a = 2, b = 1). 

III. THE GAUSSIAN MATTER FIXED POINT 

A. Minimal coupling is self consistent 

We assume that V and F are real analytic so that they can be Taylor expanded around 4> 2 = 0. A given V and F 
define a FP if the corresponding dimcnsionless potentials satisfy dtV — and dtF = 0. Because of analyticity, this 
is equivalent to requiring that all the derivatives of dtV and dtF with respect to (ft, evaluated at (ft = are zero. 
Taking n derivatives of eq. fit?]) and eq. ([2U|) with respect to 4> 2 we get 

= (d t V) in) (0) = ((d-2)n-d)V^(0) + (k~ d d t V) in) (0) ; (22) 
= (3 t F) (Tl) (0) = (n-l)(d-2)^"}(0) + (fc- (d - 2) c^) (n) (0) . (23) 

where in the last two terms the expressions in brackets can be thought of as functions of (ft 2 . We can rewrite them as 
-(k- d d t V) = k^- 2 ^- d -^—(d t V) ; -^—{k^ d ~ 2 )d t F) = k^-^-^^—idtF) . 



d{fty K ' d{<ft 2 y y ' ' d(0 2 )" W 

We now make the following Ansatz: 

V = k d X , F = k d ~ 2 £ Q , (24) 

where Ao and £o ar c numbers to be determined. This corresponds to putting a = b = in (|21|) . or in other words to 
seeting to zero all scalar self couplings. We are assuming here that all the derivatives of V and F at (ft 2 = vanish, 
so that V and F are just constants. If a FP of this type exists, we call it a Gaussian Matter Fixed Point (GMFP). 
In order to check that this ansatz defines a FP we need to show that eq. (|22[) and (f2"3"| are identically satisfied for all 

n > 1, while for n = they determine the numbers Ao and £o- 

For n > 1 the first term on the r.h.s. of eq. (|2"2")l and (|23|) vanishes because of the ansatz. There remains to show 



that gAryi (dtV) and gnpyi (dtF) are zero at (ft 2 = 0. In d = 4 one can check this explicitly by inspecting eq. (flT]) 

and eq. ()18|) . The crucial point to observe is that in dtV and dtF, whenever (ft 2 appears explicitly, it is multiplied 
by some derivative of V or F. So when the derivative removes (ft 2 , what remains is zero because of the ansatz, and 
otherwise it is zero because there remains some positive power of (ft 2 . 

In other dimensions this crucial property remains valid. In other dimensions this crucial property remains valid, 
because it is true either for the second variations (in the case of the transverse traceless tensor and transverse vector 
components) or for the matrix trace of the second variations, in the case of the scalars. Since the beta functionals 
are obtained by taking functional traces of these expressions, this property will go through for them as well i.e. for 
n > 1 the eq. (|22|) and (|23|) are identically satisfied. For a detailed proof see [ill ]. 

Thus in any dimension the ansatz works for all n > 1. There remains to solve the equations for the constant terms 
in V and F, which arc given by Ao and £q- We are going to do this numerically in section IV. In the meanwhile we 
assume that such a solution exists, and we study the properties of the linearized flow around it. 
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Linearized Flow around GMFP 



To study the linearized flow around GMFP it will be convenient to Taylor expand V and F as follows: 

oc oo 

V(ft) = J2 A 2 „ </> 2 " ; F(0 2 ) = £ ^ • (25) 



n=0 



n=0 



We define dimensionless couplings A 2 „ = fc 2 - ) "A2„ and £2,1 = k ^ 2 - )( - 1 n '^2 n , in such a way that the dimen- 

sionless potentials can be expanded as: 



f(0 2 ) = ]Ta 2 „, 



Fl 



(26) 



n=0 



n=0 



To obtain the running of dimensionless couplings we take derivatives of eq. (|19p and eq. (|2T)|) with respect to </> 2 and 
use eq. ([T(!t 



dfXvn — — 



1 <5 n au 



2\n 



r 1 S n d t F 



62=0 



n! <J(</,2)n 



(27) 



2 =O 



Because of the presence of t-derivative on the RHS of FRGE, we do not obtain the beta functions of dimensionless 
couplings directly, rather we get algebraic equations for them, solving which one get the full beta functions. 

Having defined the dimensionless couplings, we now define the stability matrix to be the matrix of derivatives of the 
dimensionless beta functions with respect to the dimensionless couplings at the FP. By definition it is a tensor quantity 
in the theory space. It will be convenient to write Vq = V and V\ = —F. One can then define the corresponding 
dimensionless potentials as V a = k d ~ 2a V a , where a is cither or 1. Then the stability matrix is given by, 



(28) 



FP 



From the above definition of the stability matrix we note that the couplings get arranged in the following order: Ao, 
Co i A 2 , £2- A4, £4 . . .. Then the matrix M at the GMFP has the following form: 



/M 00 M01 
Mn M12 
M 22 M 23 
M33 

V 



(29) 



where each entry is a 2 x 2 matrix of the form 



Mi. 



(30) 



Moreover the various non zero entries of M are related to each other by the following recursion relations (in d- 
dimensions): 



M u = (d-2)i + M ( 



where 



M 00 = 



-d 






-(d-2) 



on 



A0A0 



M 



(i + 1) (2i + 1) M01 



<5M X * 
<5M? ? 

soco 



M, 



01 



JM, x 
<5M,- x 



<5Mt ? 
5M? ? 

«0«2 



(31) 



(32) 



We can prove these facts for the one loop beta functions, i.e. neglecting the ^-derivatives of the couplings on the 
r.h.s. of FRGE. Using this we note that the running of dimensionless potentials can be written as follows: 



d t V a 



-{d - 2a) V a + (d- 2) ft V' a + H a ( V b , ft V b X, 2V b ' + Aft vf) . 



(33) 
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We have indicated that the one loop beta functional depends on <f> 2 only through the three types of combinations 
indicated as the arguments for H a . This can be verified in d = 4 by inspection of eq. (fT7|) and eq. (fl~8"|) , when one 
drops the terms proportional to dtF and dtF' in the r.h.s. The properties of the stability matrix given above follow 
by taking successive derivatives of dtV a with respect to (j) 2 at (f> 2 = 0. 

The i = entries of eq. (|28[) can be calculated by setting (f> 2 = in eq. (f33|) : 



Since 9^(0) depends only on V^(0) and V^'(0), in eq 
Mqi are given by, 



d t V a (0) = -{d - 2a) V a (0) + H a {V a (0), 2F o '(0)) . (34) 

for i = 0, only j = 0,1 will be non zero. Thus M 00 and 



(M 00 ) ab = -{d-2a)6 ab + 



5H a (V c (0), 2K'(0)) 



SV b (0) 



(Moi) ab 



SH a (V c (0), 2K'(0)) 



GMFP 

Now we take first derivative of dtV a with respect to <p 2 . This gives, 



6V<(0) 



(35) 



GMFP 



dtVl 



(d - 2a)v; + (d- 2)K' + (d - 2)~4 2 v': + + XI ^~;~A VM + 4> 2 v"Vd + 2 TC) 

dVc 6{(t) 2 V^V d ) 



SH a 



SH a 



5{2Vl + 4(f> 2 V^ 



■ (2v;" + 4y c " 



'X") 



(36) 



When we set </> 2 = 0, we note from the above equation that dtV^(0) depends only on V a (0), V^(0) and V^"(0). We 
use this in eq. (j2"5)) to calculate the i = 1 entries of the stability matrix. We note that My = for all j > 3. Now 
we find the remaining possible non zero entries. For j = 0, we note that the dependence on V a (0) is present only 

But each of these terms arc multiplied cither with V' a or V£ , 



62 =0 K<t> 2 V'V>) 



4,2=0 «(2V C '+4^V C '') 



so when we calculate the stability matrix, these terms will not contribute due to GMFP conditions (Va = for all 
i > 1). Thus we conclude that Mio = 0. 

For j = 1, we take the derivative of d t V^{Q) with respect to V b . Thus using the condition of GMFP and eq. 
we find, 



-(d-2a)6ab + (d-2)5 ab + 



SV b 



(d-2)S ab + (M 00 ) ab 



4> 2 =0 



while for j = 2 we take derivatve of <9tV^'(0) with respect to V b " /2 and use eq. (j3"5|) . Thus we get 



(^12) afc = 2 



Thus we see that for i = 1 we have, 

Mio = ; Mu = (d - 2) • 1 + M 00 ; 



•6<5 bc = 6(M 01 ) Qb 



Mi 



2-3M 



My = ,Vj >3 



(37) 



(38) 



(39) 



In order to understand the structure of the lines i > 2 we will proceed by induction. We assume that the i-th 
derivative has the following structure, 



(StVaP = -(d - 2a)V« + (d - 2) (j> 2 V^+V + i (K*)) + 



+ + 
sv c 



SH a 



8{4> 2 V'V>) 



(0 



5Ha 



(5(2 v c ' + 4^ 2 v c ") 



— - f2^ 



+ 4^ 2 v;" 



(0 



(40) 



where the (• • • ) denote expressions having at least two factors of derivatives of potentials, which are irrelevant when 
calculating the entries of stability matrix. Clearly this property is true for i = 1 . We show that if it holds for a given 



9 



value of i, then it also holds for i + Thus we take one more derivative eq. (|40]) and we find 

(d t V a )^ = -(d 2a)V^ + (d - 2) (> V^ + (i + 1) V^A +{••• + + Sl^VPrt 

K 'I oVc oVcoVd 

+ = S2 . H t - f 2 ' + f g ° - - V® f 2 V>' + 40 2 V>") ' 



<5V e <5(0 2 1/' V' V <5V r e <5(2V r c / + 40 2 y c ") 

<^ 2 K%') V ' W(0 2 u c 'uj) v v 

fr _ _ _ / . _ _ s (i) 

n a / 72 iri -Trist / 12 fW fW \ v ' 



i — (0 2 K' n')' f ^ 2 y' v3 

rti2V'VV e i> v c d 



&(<l?ViVf)8{<pVM) 

+ — = _° f 1 * ( 2 + 4 2 v;")' f 2 v c ' uj) 

<5(2 VI + 40 2 V e ") <5(0 2 VI Vfi V V 

+ - ^ - (2 V* + 40 2 V>") + ■ ■ % (2 V;' + 40 2 ^ w 

(5(2 V c ' + 40 2 VI') ^ c / <5K<*(2V r c ' + 40 2 V r c ") e V 



J(2 U e ' + 40 2 U e ")(5(0 2 U c 'U d ') 



F ^(2V> e ' + 40 2 K") W (tf 2 ^)' 

(2^ + 40 2 K'') W (2V d ' + 40 2 ^')'l • (41) 



(5(2 V c ' + 402 VI') (5(2 + 40 2 V d " 

Aside from the new terms containing two factors of derivatives of the potentials, which can be neglected for our 
purposes, the remaining terms have the same structure as eq. I|40p . Thus by induction eq. (|40| holds for all i. 
We can now use this result to calculate the entries of the stability matrix in the i-th row. Using 

(2 V' a + 40 2 VZ\ = 2 (2 i + 1) U a (m) + 40 2 V^ +2) , (42) 

we note that at = and using the condition of GMFP for calculating the stability matrix we have, 

Ma = ,Vi < (i-1) ; M u = (d-2)i+M 00 ; M M+1 = (2i+l)M 01 ; M tJ = ,Vj > (»+2) . (43) 

This completes the proof of our statements in the one loop approximation. It is difficult to extend this proof to the 
exact equation, but we see in finite truncations that the previous properties of the stability matrix remain true. 

Having established the properties of stability matrix we would like to compute its eigenvalues. The good feature 
of the block structure of stability matrix indicated in eq. (|29p is that the eigenvalues are given just by the diagonal 
blocks. Since the consecutive diagonal blocks just differ by d — 2, the eigenvalues of the consecutive diagonal blocks 
of M also differ by d — 2. This is a very strong result, because it implies that, at a GMFP, the eigenvalues of M are 
all determined by the eigenvalues of Moo- Furthermore, the off diagonal blocks of M are all determined by M01, so 
knowing Moo and M01 one can also determine all the eigenvectors. This is useful to understand the mixing among 
various operators at the FP. The smallest truncation that is required to calculate both Moo and Mqi is when we retain 
terms up to 2 in each potential. 

IV. NUMERICAL RESULTS 
A. The GMFP in d = 4. 

We now look for GMFP in various dimensions and calculate the critical exponents of the system, which are defined 
to be the opposites of the eigenvalues of M, i.e. 6j = — A», where A; is the eigenvalue. As explained in the previous 
section, it is enough to calculate the eigenvalues of M o- We do this task first in d = 4. 

In d = 4 for De-donder gauge we get the following FP equation, 

24 A ( 16tt 2 + . 1 - ) = 19, (44) 
V A - Co / 

(265 - 2304 7r 2 £oKo + (127 + 2304 7T 2 f )Ao = 3 6 ° A g . (45) 

Ao — £0 



10 



On solving these, the only real solution that we get is 

Aq = 0.00862 ; g = 0.02375 



(46) 



We now compute the critical exponents 6 of the stability matrix in this gauge. The relations given in eq. (|31j) between 
the various nonzero entries of the stability matrix are independent of the gauge. However the entires of Moo and Moi 
are gauge dependent. For d = 4 cq. (|3"Tj) reduces to, 



Ma = 2i + M 00 ; Mi.i+1 = (* + 1) (2 i + 1) M 01 
In De-Donder gauge for d = 4, the entries of Mqo are 



M l I 

Aq Aq 



32tt 2 6 



M 



A °«° 32tt 2 6 



^ {9667 + 3456tt 2 |o(169 + 2304tt 2 |o)} + A 2 , |3279 + 1152^(275 + 69127r 2 £ )} 
18A |o {551 + 128tt 2 |o(331 + 6912tt 2 |o)} 

48384A 4 , - 4435207r 2 £ 4 + A | 2 |— 9667 + 2304^(161 - 3456tt 2 | )} 



1 1 



- 3Afj |l093 + 4608tt 2 ^o(1 + 576tt 2 ^ )} + 18A§£ {551 + 1927r 2 <f (-l + 4608tt 2 £ )} 



M, 



_12 



Mr r = — 



g |-252A |o(3 + 128tt 2 | ) + 3A 2 (41 + 1536tt 2 |o) + |g(593 + 27648tt 2 £o)} 
-135A 4 , + 1309|^ + 12Ag| (145 + 1536tt 2 | ) + 36A |o(77 + 3072tt 2 | ) 



-6A^ 2 (811 + 235047t 2 | ) 
While the entries of Mqi are 



1 1 

A0A2 ~i^e 



Co(37 - 1152^o) + 2AoCo(-5 + 11527r%) - X 2 (7 + 1152tt%) 



Mx ~ p = 



3 (2A -5£ )(A -£o) 



Mi 1 

SO A2 



M, 



where 



4tt 2 9 
24f (Ao - |o) 2 

e 

72g (A - go) 2 



9 = -18Ao6(l + 128tt%) + 3Ag(5 + 384tt 2 £ ) + ^(-17 + 11527r 2 £ ) 



(47) 



(48) 



(49) 
(50) 



(51) 

(52) 

(53) 
(54) 
(55) 

(56) 



The relations eq. (|4"T)) tells that the critical exponents of consecutive diagonal blocks will differ by 2. In the 
truncation where we keep terms till </> 2 in each potential, the critical exponents are, 

2.143 ±2.879i , 0.143 ± 2.879i (57) 

The critical exponents 2.143 ± 2.879i correspond to eigenvalues of Moo, while the critical exponents 0.143 ± 2.879z 
which are shifted by 2 correspond to the eigenvalues of M\\. This justifies our claim. The eigenvectors in this 
truncation are 

0.3557 ± 0.3776i \ / (-18.059 ± 7.31(H) x 10~ 4 

0.8549 I I (-30.723 ± 10. 763i) x 10~ 4 

' 0.3557 ±0.3776i 

/ V 0.8549 



(58) 
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FIG. 1: The gauge dependence of GMFP and critical exponents in d — 4. For convenience of depicting the values in plot, we 
have rescaled Aq and £o by 1000 and 100 respectively. 



where the first complex conjugate pair of eigenvector correspond to critical exponents 2.143 ± 2.879i, while the second 
pair correspond to critical exponents 0.143 ± 2.879i. 

We then looked for GMFP in other gauges. In d = 4 we consider various values of the gauge parameters a and 
p. To study the gauge dependence we considered 50 different values of a in the range to 1.225 at step of 0.025, 
and 25 different values of (3 in the range —1 to 1.4 at interval of 0.1. For each combination of a and (3 we solved the 
FP equation obtained for Ao and £o- I n general, this produces a set of FPs. In order to choose the correct GMFP 
from that set, we plot all the real FPs to see which one is continously followed in other gauge values and which ones 
are spurious. For example one can take any value of /3, and plot all the real FPs for various values of a. Some FPs 
don't exits for all values of a, and are assumed to be truncation artifacts. Only one GMFP exists for all values, and 
is continuous. This observation of continuity in a and j3 is useful to write a code for selecting the right GMFP for 
various gauges. After calculating the GMFP we calculate the critical exponents of Moo- We then plot the GMFP and 
critical exponents against the various gauge values and generate 3D graphs. In d = 4 we obtain the graphs shown in 
Fig. (JT]). We note that the existence of the FP has been actually verified in a much larger range of values of a and (3. 



B. The GMFP in other Dimensions 



We now look for the GMFP in other dimensions. For any d > 2, in Dc-donder gauge, the FP equation for Ao and 
fo is given by, 

2dX + (3d-2)Ao + ( d -2)(^ + d-l)| = 
(47r)^r(2 + f)(A -£o) 



A\ 2 + B\qZ + C% 
(Ao - fo) 



+ D(A -£o) = , (60) 



where 



A = -(d-l)(d 3 + 2d 2 + 36d + 24) , 

B = -(d-l)(d 5 -17d 3 -38d 2 -96d-48) , 

C = (d 6 - 13d 5 + 32d 4 - 104d 3 + 72d 2 + 36d + 24) , 

D = 24d(d-l)(d-2)(4^) d/2 r^2 + . (61) 

Solving these equations we find that in other dimensions, it is possible to have more than one real solution. But when 
we plot all the real the solutions against various d in a graph, we notice that not all solutions exist in all dimensions. 
Only one solution exists in all dimensions, and is continuous in d. Besides, the ones which don't exist in all dimensions, 
have large critical exponents and are probably unphysical. In Fig. 2 we plot the position of the GMFP for 2 < d < 11, 
both in terms of Ao and £o an d of the more familiar dimcnsionlcss cosmological constant and Newton constant 

Ao = ; fo = ■ (62) 

16ttG 16ttG 

After having found the GMFP in various dimensions, we set to calculate their critical exponents. In arbitrary 
dimensions, the various blocks of the stability matrix obey cq. (f3"T|) . We plot the critical exponents of Mqq for various 
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FIG. 2: In the first graph we plot the GMFP Aq and £ { * in various dimensions. In the second plot we calculate the corresponding 
FP values of the cosmological constant A* and Newton's constant G* in various dimensions. 




2 4 6 8 10 2 4 6 S 10 



d d 

FIG. 3: Critical exponents at the GMFP in various dimensions. The left panel shows the real part of the critical exponents, 
the right panel shows the imaginary part of the critical exponents. We note that below d = 2.8 the critical exponents becomes 
real. 

dimensions. From the graph Fig. ([3]) we note that around d = 2.8 there is bifurcation. Below d < 2.8 the critical 
exponents are no more complex. 

A summary of the properties of the GMFP in various dimensions is given in table (JTJ) . Notice that for all the 
dimensions considered, the real part of the critical exponents is greater than d — 2 and less than 2(d — 2). As a result, 
in all these cases there are exactly two pairs of complex conjugate critical exponents with positive real part, i.e. four 
relevant directions. 



d 


AS 




A* 


G* 


0i 


02 


2.001 


4.968 


xlO" 


b 


2.386 


xlO 


i 


1.041 


xlO" 


. 


8.339 xl0 _b 


2.001 


0.001 


3 


1.605 


xlO" 


2 


1.047 


xlO" 


-l 


7.666 


xlO" 


2 


1.900 xlO" 1 


1.627 + 0.754 i 


1.627 


- 0.754 i 


4 


8.620 


xlO" 


3 


2.375 


xlO" 


-2 


1.814 


xlO" 


1 


8.375 xlO" 1 


2.143 + 2.879 i 


2.143 


- 2.879 i 


5 


2.669 


xlO" 


3 


5.744 


xlO" 


-3 


2.323 


xlO" 


1 


3.463 


3.236 + 4.996 i 


3.236 


- 4.996 i 


6 


6.230 


xlO" 


4 


1.207 


xlO" 


-3 


2.581 


xlO" 


1 


1.648 xlO 


4.818 + 7.039 i 


4.818 


- 7.039 i 


7 


1.225 


xlO" 


4 


2.235 


xlO" 


-4 


2.740 


xlO" 


1 


8.900 xlO 


6.744 + 9.004 i 


6.744 


- 9.004 i 


8 


2.133 


xlO" 


5 


3.738 


xlO" 


-5 


2.853 


xlO" 


1 


5.322 xlO 2 


8.945 + 10.904 i 


8.945 


- 10.904 i 


9 


3.380 


xlO" 


6 


5.747 


xlO" 


-6 


2.941 


xlO" 


r 


3.462 xlO 3 


11.396 + 12.748 i 


11.396 


- 12.748 i 


10 


4.960 


xlO" 


7 


8.228 


xlO" 


-7 


3.014 


xlO" 


l 


2.418 xlO 4 


14.089 + 14.537 i 


14.089 


- 14.537 i 


11 


6.817 


xlO" 


8 


1.107 


xlO" 


-7 


3.079 


xlO" 


r 


1.797 xlO 5 


17.025 + 16.261 i 


17.025 


- 16.261 i 



TABLE I: Position of GMFP and critical exponents for various dimensions. 
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FIG. 4: The gauge dependence of GMFP and critical exponents in d — 6. For convenience of depicting the values in plot, we 
have rescaled Ag and £q by 10 4 and 10 3 respectively. 

Finally we studied the gauge dependence in different dimensions in the same wasy as we did in d = 4, for example 
in d = 6 we obtain the graphs shown in Fig. ([¥]). 

V. OTHER NON TRIVIAL FIXED POINTS 

Having discussed the existence and properties of the GMFP, we can ask ourselves whether there exist other FP 
where the scalar field has nontrivial self-interactions. We look for (truncated) polynomial FP potentials of the form 

a b 

V(0 2 ) = J2 4> 2n ; F(^) = J2 ^ ■ ( 63 ) 

with finite a > 1, b > 0. Such potentials are known not to exist in a pure scalar theory in four dimensions [431 ] . so we 
consider it unlikely that they exist in the presence of gravity, In fact the outcome of our numerical searches is that 
no such FP's appear to exist in dimensions 4, 5 and 6. (Some FP do appear in certain truncations but not in others, 
so they are likely to be just truncation artifacts.) 

The situation is somewhat different in three dimensions. We know that pure scalar theory in d = 3 has the Wilson- 
Fisher FP [44|]. This FP can be seen in our calculations by taking the limit G — > (where Newton's constant G is 
related to £0 = l/167rG) and Ao — ► 0, in which case gravity decouples. Solving the FP equations of the scalar field 
in the LPA, truncated to order 4 , one gets A2 = —0.0385 and A| = 0.3234, with critical exponents 9\ = 1.843 and 
62 = —1.176. (These are not very good values, but we quote them here for the sake of comparison with what we find 
in the presence of gravity.) The FP persists when one goes to higher truncations. 

One wonders whether there exists a "gravitationally dressed" Wilson-Fisher FP, with nonvanishing G, namely a 
FP where gravity and the scalar simultaneously have nontrivial interactions. Again in certain truncations one finds 
various FPs which turn out to be truncation artifacts. There seems however to exist one genuine FP: we find it in 
all truncations where a > b, and it has very similar properties in all truncations. To explore its properties we have 
looked in two directions: increasing simultaneously a and b, or keeping 6 = and increasing a. 



(a,b) 


AS 


A^ 


A2 


A^ 


As 


Co* 


C2* 


C2 




Cs* 


(2,1) 


0.0196 


-0.1646 


-0.1595 






0.1088 


-0.03108 








(3d) 


0.01994 


-0.1758 


-0.1958 


-0.2796 




0.1096 


-0.03810 








(4,1) 


0.02002 


-0.1783 


-0.2041 


-0.3466 


-0.5579 


0.1098 


-0.03969 








(2,2) 


0.01894 


-0.1408 


-0.1241 






0.1071 


-0.01122 


0.04297 






(3,2) 


0.01971 


-0.1680 


-0.1848 


-0.2879 




0.1089 


-0.03131 


0.01731 






(4,2) 


0.01988 


-0.1735 


-0.1975 


-0.3544 


-0.5687 


0.1093 


-0.03542 


0.01121 






(3,3) 


0.01911 


-0.1469 


-0.1469 


-0.1935 




0.1074 


-0.01420 


0.05017 


0.1617 




(4,3) 


0.01953 


-0.1618 


-0.1768 


-0.3083 


-0.6569 


0.1084 


-0.02571 


0.03197 


0.1102 




(4,4) 


0.01923 


-0.1512 


-0.1572 


-0.2496 


-0.4911 


0.1077 


-0.01728 


0.04732 


0.1765 


0.3868 



TABLE II: Position of Nontrivial FP in d — 3 for various truncations. 

Tables ((IT]) and (|III[) give the position and critical exponents of this FP for a > b and b < 4. One notices that 
A*!ri < f° r all n > in the table. It is computationally demanding to continue in this direction, so to have some 
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(a, b) 


o[ 


0i 


03 


04 


05 


06 


07 








(2,1) 


1.648 


0.592 


-0.956 


-3.902 


-13.46 












(3,1) 


1.650 


0.554 


-1.079 


-3.776 


-11.20 


-29.397 










(4,1) 


1.650 


0.543 


-1.105 


-3.673 


-10.02 


-24.01 


-49.31 








(a, b) 


0'i 


0i 


02 


02 


05 


0c 


07 


08 






(2,2) 


1.649 


0.656 


-7.979 


1.261 


-0.559 


-3.192 










(3,2) 


1.652 


0.589 


-7.933 


3.909 


-0.835 


-3.578 


-27.67 








(4,2) 


1.652 


0.570 


-7.635 


4.083 


-0.898 


-3.626 


-22.64 


- 47.78 






(a, b) 


777 


7T77 

0i 


777 

02 


7T77 

02 


777 

03 


7777 

03 


07 


08 


09 




(3,3) 


1.649 


0.641 


-6.703 


2.097 


-14.12 


8.990 


-0.512 


-2.991 






(4,3) 


1.651 


0.603 


-6.448 


3.343 


-13.94 


10.91 


-0.657 


-3.287 


-42.28 




(a, b) 


05 


0i' 


02 


2 ' 


03 


<% 


03 


<% 


09 


010 


(4,4) 


1.650 


0.630 


-5.958 


2.008 


-12.88 


7.966 


-20.07 


19.03 


-0.513 


-2.977 



TABLE III: Critical exponents at Non trivial FP in d = 3 for various truncations. When the critical exponents are complex 
we write them in the form 9' e ± i8'/ 

indication on the sign of A^ n for higher n we considered a simple truncation where F is constant, i.e. £ n = for 
?i > 0. In this case we could push the truncation up to a = 8. The results are given in tables (|IV[) and ((Vj. One sees 
that the coefficients of the potential are indeed all negative. Furthermore, the coefficients grow in absolute value, so 



(a,b) 


AS 


A 2 


A3 


A 6 


A« 


Aio 


At 2 


Al4 


Aio 


Co* 


(1,0) 


0.01813 


-0.1088 
















0.1060 


(2,0) 


0.01880 


-0.1343 


-0.1561 














0.1065 


(3,0) 


0.01894 


-0.1395 


-0.1942 


-0.2633 












0.1066 


(4,0) 


0.01898 


-0.1407 


-0.2032 


-0.3284 


-0.4998 










0.1066 


(5,0) 


0.01899 


-0.1410 


-0.2053 


-0.3437 


-0.6182 


-0.9604 








0.1066 


(6,0) 


0.01899 


-0.1411 


-0.2058 


-0.3472 


-0.6452 


-1.180 


-1.826 






0.1066 


(7,0) 


0.01899 


-0.1411 


-0.2059 


-0.3479 


-0.6511 


-1.228 


-2.229 


-3.380 




0.1066 


(8,0) 


0.01899 


-0.1411 


-0.2059 


-0.3481 


-0.6524 


-1.238 


-2.313 


-4.091 


-5.977 


0.1066 



TABLE IV: Position of Nontrivial FP in d = 3 for other truncations. 



(a,b) 


0i 


91 


03 


04 


05 


06 


07 


08 


09 


010 


(1,0) 


1.659 


0.753 


-1.699 
















(2,0) 


1.675 


0.745 


-1.594 


-10.59 














(3,0) 


1.679 


0.742 


-1.485 


-8.384 


-22.41 












(4,0) 


1.68 


0.741 


-1.434 


-7.341 


-18.19 


-36.96 










(5,0) 


1.68 


0.741 


-1.414 


-6.84 


-15.93 


-31.12 


-54.10 








(6,0) 


1.68 


0.741 


-1.407 


-6.609 


-14.70 


-27.47 


-47.25 


-73.60 






(7,0) 


1.68 


0.740 


-1.405 


-6.509 


-14.02 


-25.287 


-42.062 


-66.663 


-95.172 




(8,0) 


1.68 


0.740 


-1.405 


-6.469 


-13.67 


-23.94 


-38.77 


-59.72 


-89.58 


-118.4 



TABLE V: Critical Exponents at Non trivial FP in d = 3 for other truncations. 

the series for V has a very small radius of convergence. This is similar to the situation discussed in (43j . making the 
FP unphysical. So we conclude that also in three dimensions there is probably no physically viable FP besides the 
GMFP. 



VI. CONCLUSIONS 

The results of this paper confirm and extend the findings of (35|. The GMFP is found to exist also in other 
dimensions and in other gauges, and (with the possible exception of d = 3) there does not seem to be other FP's with 
nontrivial scalar self-interactions. In four dimensions this agrees with the findings of [43}. 

These results may be applied in various settings. The beta functions given in Appendix A contain the full dependence 
on the dimensionless parameters Aq, £o, A2, £2, A4, without making any assumption on the value of these couplings 
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(which in the case of the first three means the ratio between the dimcnsionful couplings Ao, £o, A 2 and the RG scale 
k). In particular, threshold effects are taken into account by the denominators 1 + 2A2 and £0 — Ao- One can easily 
recognize among various terms the ones that are obtained in perturbative approximations, but we emphasize that the 
derivation of these beta functions using the FRGE does not require that the couplings be small. 

The most natural application of these results seems to be in the context of early cosmology, where a scalar field is 
used to drive inflation. In an asymptotic safety context, it would be attractive to obtain inflation as a result of FP 
behavior along the lines of [23| . In fact the energy scale involved is sufficiently high that one may expect quantum 
gravity effects to play some role. Alternatively, it would also be of interest to apply the flow equations derived here 
to the scalar tensor theory, e.g. to improve the results of [37| . 

According to various speculations, quantum effects may play a role also on very large scales, and then again the 
RG flow of scalar-tensor theory could become relevant. In this connection we recall that scalar-tensor theories of a 
different type also arise in the conformal reduction of pure gravity, and have been studied from the FRGE point of 
view in (Tall]. 

We have mentioned in the Introduction that scalar-tensor theories can be reformulated classically also as pure 
gravity theories with f(R) type actions, and one may wonder whether there is a relation also between their RG flows. 
In particular one could ask whether the FP that was found in [l3l [l4| has a counterpart in the equivalent scalar-tensor 
theory. At first sight one would think that this is not the case, because the choice of cutoff breaks the classical 
equivalence between these theories. Still, this point deserves a more detailed investigation. 

Another direction for research is the inclusion of other matter fields. As discussed in the introduction, if asymptotic 
safety is indeed the answer to the UV issues of quantum field theory, then it will not be enough to establish asymptotic 
safety of gravity: one will have to establish asymptotic safety for a theory including gravity as well as all the fields that 
occur in the standard model, and perhaps even other ones that have not yet been discovered. Ideall y o ne would like to 
have a unified theory of all interactions including gravity, perhaps a GraviGUT along the lines of 45]. More humbly 
one could start by studying the effect of gravity on the interactions of the standard model or GUTs. Fortunately, 
for some important parts of the standard model it is already known that an UV Gaussian FP exists, so the question 
is whether the coupling to gravity, or some other mechanism, can cure the bad behavior of QED and of the Higgs 
sector. That this might happen had been speculated long ago [HI; see also [46| for some detailed calculations. It 
seem that the existence of a GMFP for all matter interactions would be the simplest solution to this issue. In this 
picture of asymptotic safety, gravity would be the only effective interaction at sufficiently high scale. The possibility of 
asymptotic safety in a nonlinearly realized scalar sector has been discussed in [471 ]. Aside from scalar tensor theories, 
the effect of gravity has been studied in fill for gauge couplings and [4i| for Yukawa couplings. 



APPENDIX A: EXPLICIT BETA FUNCTIONS 



In section II. B we have presented equations which in principle determine the beta functionals for V and F in d = 4 
and in the gauge a = and (3—1. At one loop the beta functionals can be immediately read off from there, but if one 
wants to get the "improved" beta functionals (meaning that no approximation is made beyond the truncation), then 
the system is too complicated to be solved. It can be solved if we assume that V and F are finite polynomials. We 
give here a set of linear equations that determine the beta functions in a five coupling truncation, including Ao, £oj A2, 
£21 A4. These equations can be easily solved using algebraic manipulation software. This exercise will also enable us 
to compare with familiar one loop results. These beta functions had been written previously in [H|, but since there 
we had left the cutoff generic, it was not possible to compute the integrals over momenta, which are contained in the 
expressions Q2 and Q\. 1 Here the integrals have already been performed, using an optimized cutoff [36j, so the beta 

functions are in closed form and completely explicit. We use the notation rj = = + 2. 



<9 t A 
dtto 



4An 



2£o 



1 



32tt 2 
1 

384tt 2 



1 



1 + 2A 2 
4 



6An 



£0 



Ao 
24£ 2 



V 5£o-2A 
96^ fo-Ao ' 

, 8|o(7fo-2A ) 



1 + 2A 2 (1 + 2A 2 ) 2 



(€0 - Ao) 2 



rj 17£ 2 + f 8£ A - 1 5A 2 



11527T 2 



(£0 - Ao) 5 



(Al) 
(A2) 



1 Note that the notation used in [35l | is opposite to the one used here: parameters with a tilde are dimensionful, those without tilde are 
dimensionless. The beta functions written in l35t contain a number of transcription errors, which however do not affect the subsequent 
results. 
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If one neglects the terms involving r\ and 3t6 m t ne r.h.s., then the remaining terms are the one loop beta functions 
for the couplings. One can recognize among them some familiar terms. The term containing — I8A4 in the first line 
of eq. (|A3[) and the term containing 216A 2 in the third line of eq. (|A5[) are the familiar beta functions of the mass 
and of the coupling in </> 4 theory in flat space. The term containing 72A 4 (1 + 126) m t ne third line of eq. (|A4[) is also 
known from earlier calculations [37l [HI . 

Notice the ubiquitous appearance of the factors 1/(1 + 2A 2 ), which represent threshold effects for the contributions 
of scalar loops: for k 2 ^> A 2 , A 2 -C 1 and the denominator can be approximated by 1, whereas for k 2 -C A 2 , A 2 ^> 1 and 
the term is suppressed. The denominators 6 — Ao have a somewhat similar effect. When written in terms of the more 
familiar variables A and G defined in eq. (f62|) . they give rise to denominators (1 — 2A). These can be approximated 
by 1 when A k 2 but they vanish when the dimensionless cosmological constant A tends to 1/2, corresponding to 
an infrared singularity in the RG trajectories. This is well documented in the literature. 

The term proportional to 96A4 in the third line of eq. (|A5[) is the leading gravitational correction (of order G) 
to the running of the scalar self coupling. Note that for small Ao, A 2 and 61 the denominator and the bracket to 
its right can be expanded as 1 + O(Ao) + 0(A 2 ) + 0(6)- 2 The order of magnitude and sign of this term agree 
with the calculations done in (34|, in the gauge a = 0. One should not expect the results to agree exactly, because 
this term is gauge dependent and the calculation was done in a different gauge (namely f3 = —1). Notice that this 



2 on a related note, we also observe that the first term in dt\4, which is proportional to G 2 , vanishes when we set A2 = and £2 = 0. 
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term is proportional to A4, thus when we set A4 to zero, the beta function for A4 does not get any contribution from 
gravity, in agreement with the general statement that minimal coupling is self consistent. We observe that the same 
phenomenon happens in the case of the Yukawa coupling [4!| and of the gauge coupling [48| . 
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